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INTRODUCTION 

The  dynamics  of  crystal  lattices  have  been  of  interest  for  some  time, 
especially  with  regard  to  the  vibrational  motion  of  the  lattice  points.  In 
1912  Debye  (l)  obtained  a  frequency  distribution  by  treating  the  solid  as  a 
continuum.  At  that  time  the  specific  heat  obtained  from  the  Debye  frequency 
distribution  agreed  veil  with  experiment;  however,  the  theory  was  in  disagree- 
ment with  the  known  periodic  structure  of  crystals.  Born  and  yon  Karman  (2) 
in  1912  examined  the  dynamics  of  a  lattice  of  discrete  particles,  however 
their  method  was  computationally  more  complicated  than  the  Debye  theory  and 
the  simpler  theory  continued  in  use  until  its  inability  to  fit  experimental 
data  at  low  temperatures  became  apparent.  Blackman  in  1935  (5)  and  1937  (^) 
revived  the  Born- von- Karman  treatment,  and  his  work  indicated  that  the 
Born-von-Karmaja  frequency  distribution  had  the  characteristics  necessary  for 
better  agreement  with  experiment.  Blackman'a  calcvilations  for  a  two 
dimensional  lattice  were  later  shown  to  be  correct  by  exact  calculations  of 
Montroll  (5)  and  Bowers  and  Rosenstock  (6).  The  revival  of  interest  in  the 
Born-von-Karman  method  has  stimulated  theoretical  investigations  on  the 
general  characteristics  of  the  frequency  distribution  by  such  workers  as 
Van  Hove  (7)  and  Phillips  (8)  among  others. 

With  the  return  of  the  non-continuum  treatment  of  crystal  lattices  the 
question  of  the  types  of  interactions  between  lattice  points  became  pertinent. 
Perhaps  the  most  common  method  for  cubic  crystals  assumes  central  forces 
between  nearest  and  next  nearest  neighbors.  This  central  force  model  was 
used  by  Blackman  in  his  calculations  for  two  dimensional  lattices.  It  has 
also  been  used  by  Bowers  and  Rosenstock  (6),  and  by  Montroll  (5)  in 


obtaining  exact  frequency  distributions  for  two  dimensional  lattices.  The 
central  force  model  was  used  by  Jlne  (9)  and  later  Clark  (lO)  in  obtaining 
frequency  distributions  for  body  centered  cubic  crystals,  and  by  Leighton  (ll) 
for  face  centered  cubic  crystals.  Recently  Shaniland  (12)  has  used  the 
centraJL  force  model  for  both  face  centered  and  body  centered  cubic  crystals. 
Montroll  and  Peasly  (ij)  used  central  forces  in  calculating  the  frequency 
distribution  of  body  centered  cubic  crystals  using  Montroll 's  (l^)  moment 
method. 

Although  it  has  been  customary  to  assume  a  central  force  model,  this 
assvuuptlon  leads  to  the  Cauchy  relationship  between  the  elastic  constants  of  the 
crystal,  a  relationship  which  is  not  born  out  by  experiment.  Owing  to  the  pre- 
diction of  the  Cauchy  relationship  from  the  use  of  central  forces  alone  various 
models  incorporating  non-central  forces  have  been  examined.  Smith  (15)  pro- 
posed a  model  for  the  dieimond  lattice  including  arbitrary  forces  between  nearest 
neighbors,  and  central  forces  between  next-nearest  neighbors.  Rosenstock.  and 
Newell  (16)  have  examined  a  simple  cubic  lattice  using  non-central  forces 
between  nearest  neighbors,  and  Singh  and  Bowers  (17)  have  examined  a  body 
centered  cubic  crystal  using  non- central  forces  between  nearest  neighbors  and 
central  forces  between  next-nearest  neighbors.  This  paper  treats  a  body 
centered  cubic  lattice  using  non-central  forces  of  the  type  Introduced  by 
Gazls,  Herman,  and  Wallis  (18)  In  treating  surface  waves  in  a  simple  cubic 
crystal.  In  this  model  non-central  forces  consisting  of  angular  stiffness 
forces  associated  with  primitive  angles  of  the  lattice  are  included  along 
with  the  central  forces  between  nearest  and  next-nearest  neighbors.  The 
addition  of  this  third  force  constant  makes  it  possible  to  derive  explicit 
relationships  between  the  elastic  constants  of  the  crystal  and  the  lattice 


model  force  constants.  In  addition  to  the  angles  considered  by  Gezls  et.  al., 
additional  angles,  formed  with  one  side  the  nearest  neighbor  distance  and  the 
other  side  the  next -nearest  neighbor  distance,  were  Included  In  this  work. 
The  addition  of  the  second  angular  terms  resulted  in  four  force  constants 
which  were  matched  with  the  three  elastic  constants  by  assuming  a  relationship 
between  the  two  angular  force  constants. 

As  is  usual  plane  wave  solutions  to  the  equations  of  motion  are  assvuned 
resulting  in  a  secular  equation  in  the  square  of  the  frequency.  The  roots  of 
this  equation  were  fovind  at  approximately  59,000  solution  points  within  the 
first  l/U8th  Brillouin  zone  using  an  IBM  709O  digital  computer.  Frequency 
distribution  histograms  were  constructed  by  determining  the  number  of  modes 
of  oscillation  within  specified  frequency  intervals.  When  the  non-central 
forces  are  ignored  this  model  reduces  to  the  ordinary  central  force  model, 
and  in  this  case  the  calculations  agree  with  the  previous  work  of  Fine  (9) 
and  Clark  (lO). 

One  of  the  primary  aims  of  the  Investigation  was  to  examine  the  effect  of 
the  Inclvislon  of  non-central  forces  on  the  frequency  distribution.  Comparison 
between  the  central  force  model  and  the  model  employed  here  was  made  for 
tungsten  which  is  essentially  isotropic  and  vanadium  which  is  one  of  the  more 
anisotropic  body  centered  cubic  crystals.  The  calculations  with  vanadium  were 
also  compared  with  the  frequency  distributions  obtained  by  the  slow  neutron 
experiments  of  Steward  and  Brockhous  (19),  and  Eisenhauer  et.  aL.  (20). 
Dispersion  curves  along  specified  directions  were  calculated  in  order  to 
compare  the  model  with  the  experimentally  determined  dispersion  curves  of 
Low  (21)  for  iron. 


DERIVATION  OF  UTTICE  EQUATIONS  OF  MOTION 

The  equations  of  motion  for  the  Born- von- Karman  model  of  the  crystal 
lattice  are  required  to  compute  the  vibrational  modes  of  a  crystal.  The 
nature  of  the  interactions  between  the  lattice  points  must  be  assumed. 
Perhaps  the  most  common  assvunptlon  Is  the  approximation  that  the  only  inter- 
actions are  central  forces  between  nearest  and  next-nearest  neighbors  of  the 
crystal  lattice.  The  model  used  here  for  a  body  centered  cubic  lattice  of 
identical  particles  assumes  non-central  angular  stiffness  forces  between 
primitive  and  secondary  angles  of  the  lattice  as  well  as  the  usual  central 
forces  between  nearest  and  next-nearest  neighbors. 

Figure  1  shows  the  <jonflguratlon  of  the  lattice  points  and  the  indexing 

of  these  points.  A  vector  S.  will  refer  to  a  general  infinitesimal  displace- 

J 

ment  of  the  Jth  atom  in  the  lattice.  The  displacement  vector  S.  has 

J 

components  (u,,  v  ,  w  )  in  the  x,  y,  and  z  directions.  Although  the  equations 
of  motion  for  central  forces  appear  elsewhere,  for  example  in  DeLaunay  (22), 
they  will  be  repeated  here  in  order  to  provide  a  continuity  of  notation. 
The  central  force  between  nearest  neighbors  may  be  represented  by 

where  a  is  the  nearest  neighbor  force  constant  and  7.  is  a  unit  vector  with 

J 

components  (n  ,  v  ,  tj  )  with  direction  along  the  line  of  action  between 

/  5 
the  central  atom  and  the  nth  atom  and  the  sum  is  over  all  nearest  neighbors 

(nn).  The  components  i^^,   v^,   r\^  take  on  the  values  +  1  depending  on  the 

location  of  the  nth  atom.  Expanding  the  scalar  product  in  Eq.  (l),  we  may 

write  the  force  acting  on  the  central  atom  in  the  x-direction  as  follows: 


EXPIAMTION  OF  PIATE  I 

Fig.  1.  The  indexing  of  a  body  centered  cubic  lattice 
using  1/2  the  lattice  spacing  as  the  basic 
unit. 
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EXPLANATION  OF  PLATE  II 

Fig.  2.   Body  centered  cubic  lattice  showing  primitive 
(0  )  and  secondary  (0  )  angles. 

Fig.  3.   Illustration  of  the  two  triangles  considered 
in  deriving  the  angular  stiffness  forces. 


PLATE  1 1 


l+l.m-l,n+l 
Figure  2 


l+2,m,n 


Figure  3 


72  secondary  angles  of  the  type  0^  in  Fig.  2  having  a  force  constant  /_  .  To 
find  the  potential  function  given  by  Eq.  (5)  it  is  necessary  first  to  construct 
the  vectors  n.  In  order  to  construct  these  vectors  only  the  indexing  of  the 
tvo  Identical  triangles  shown  in  Fig,  3  need  he  examined.  In  this  figure  the 
point  labeled  o  refers  to  the  center  atom,  m  and  k  to  atoms  on  the  cube  corners, 
and  q  to  a  next-nearest  neighbor  atom.  The  vectors  x,  ,  x  ,   and  x  ,  each 
drawn  from  the  origin  to  the  atom  vmder  consideration,  are  conveniently  found 
by  taking  their  components  to  be  the  differences  of  the  indices  of  the  kth, 
rath,  or  qth  atom  from  those  of  the  center  atom,  namely  ^,   m,  n.  For  example. 
If  the  kth  atom  is  the  -t+l,  m+1,  n+1  lattice  point  then  x.  has  all  components 
equal  to  one.  The  vectors  n  may  now  be  taken  as  a  Uneeir  combination  of  the 
vectors  x.  and  x  in  triangle  kom  which  will  be  called  a  type  1  triangle,  and 
-x  and  (x  -X  )  in  triangle  omz  which  is  designated  type  2.  There  is  the  added 

AU  ill     \^ 

condition  that  the  vectors  n  be  normal  to  the  sides  of  the  angles  and  in  the 
direction  of  decreasing  angle.  It  should  be  mentioned  that  the  choice  of  the 
direction  of  the  n  is  arbitrary  but  must  be  consistent.  The  derivation  of  the 
non-central  force  on  the  center  atom  will  be  carried  out  for  triangle a  of  type 
1,  only,  as  a  simple  transformation  will  give  the  contribution  dixe  to 
triangles  Including  a  next-nearest  neighbor  atom.  The  force  on  the  central 
atom  due  to  angular  deformations  of  the  first  type  triangle  may  be  written  as 

F  a  -  •=  V  V 
1    2   o  *1 

""^  Ty^l  ^^o^WV^o^-V^l^W-mo)' 
+  72  (C2(V\)-%n^l(V\)*\o^^l  ^ 
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where 

o 


^-  =  ^  ^   -^  J   57-   -^  ^  37- 


0  0  0 

In  this  equation  C^  and  Cp  are  constants  which  take  into  account  the  difference 
in  the  leiagths  of  the  sides  of  the  angles  and  the  magnitudes  of  the  normal 
vectors • 

The  normal  vectors  n.  ,  n  ,   and  il   are  defined  hy  the  equations 

so   "  m  '    K. 
\o  =  ^  ■  ^  ^m  ^'^^ 

where 

e  e  — — — —  £3  cos  0  t  (8) 

The  change  in  potential  energy,  Vp>  due  to  angular  changes  in  type  2 
triangles  may  be  written  directly  from  the  potential  function  V,  in  Eq,.  (6) 
as  the  triangles  are  identical  except  for  labeling.  Referring  to  Fig.  5, 
the  potential  function  Vp  is  obtained  from  V  by  the  following  transformations 
of  indices 

m— ^  p  , 

k— >  q  , 

P— >  a  I 
and  by  letting 

(9) 


u 


The  expansions  of  Eqs.  (6)  and  the  resulting  expression  with  the  substitutions 
given  In  (9)  give  the  ftngi1«T'  forces  in  the  x-dlrection  in  terms  of  displace- 
ments of  the  lattice  points  in  the  x,  y,  and  z  directions.  The  expansions  are 
quite  tedious  hut  straightforward.  The  total  contribution  from  the  angular 
term  results  in  a  force  on  the  central  atom,  in  the  x  direction,  given  hy: 

1  •  F.  «7t  C?  fup   „  [-52(l-e)^-8(5-2e+3e2)+8(l-e)%-8e(l-e)(V-fW) 

44(l+e)2(u^2^^^^+u^2^^^^)-m-e)2(u^^^^2,n-*-^^,m-2,n+%m,n+2+Ht,„,,n-2^^J 
"■^2  {H,m,n^-2cf(>2e+3e2)^C3_C2(l-e)-l6c|]  ^^^^ 

+U[  2cJ(3-2e+5e^ )+(V+W) [ 2C^(5e^-2e-l)+8c^C2(l-te )] 

432c|-l6c^C2(l.e)](u^^^2^^+u^^^_2,n-^%m,n+2+^^,m,n.2)3    > 
where 

^  "   ^^tfl,m+l,n+l"*'^-tfl,m+l,n-l'*''MU-l,m-l,n+l'^^-tfl,ii>-L^4>.l* 
^-t.l,m+l,  n+l'*'^-t.l,  m+1,  n-l'^^l^m-l,  n+l'^-^1,  m- 1,  a-l) 

^  ■  ^"^^l,m+l,n+l'*'^-tn,mfl,n-l"'^-tfl,m-l,n+l"^-tfl,m-l,n-l"  (ll) 

"^-Ul,  m+1,  n+l""^-^l>  m+1,  n-1+^-t-l,  m-1,  n+l'^^-t-i,  m-1,  n-l) 

*  "  ^^•tn,m+l,n+l'^'tn,m+l,n-l"*'^-tfl,m-l,n+l"^-tfl,m-l,n-l" 

^-^1,  mH,  n+l'^^-t-l,  m+1,  n^l"^-t-l,  m-1,  n+l'*^-t-l,  m-1,  n-1^  * 

For  small  relative  displacements  the  motion  of  the  other  lattice  points 
may  "be  related  to  the  displacement  of  the  central  atom  and  its  spatial 
derivatives  hy  a  Taylor  series  expansion  about  the  central  atom.  In  the  limit 
of  small  displacements,  however,  the  coefficient  of  the  first  order  term 
u^  jjj  ^  must  tend  to  zero.  This  condition  is  satisfied  by  Eq,.  (lO)  when  e  is 
replaced  by  its  value,  l/j,  as  found  from  Eq.  (6).  It  is  also  desirable  that 


12 

the  non-central  forces  discussed  here  he  Invariant  vmder  rigid  hody  rotations. 
Equation  (lO)  satisfies  this  condition  for  the  ratio  C,/Cp  =  l/2,  which  takes 
into  account  the  difference  in  the  lengths  of  the  sides  of  the  angles  and  the 
different  magnitudes  of  the  normal  vectors. 

The  equation  for  the  force  in  the  x-direction  on  the  central  atom  due  to 
the  central  and  non-central  forces  considered  here  may  now  be  written  as 


X 


m  \io 


^  't,ra+2,n  'C,m-2,n  't,m,n+2  't,m,n-2 

+  U  [f  +  2  7^^  (|§)  +  5  72  if)]   +  (V4V)  [f  -  7i  if)   +1  72  if)]   , 

where  U,  V,  emd  W  are  defined  as  in  Eq.  (ll).  The  forces  in  the  y  and  z 
directions  on  the  central  atom  are  found  hy  a  circular  permutation  of  u,  v, 
and  w. 


(12) 
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SOLUTION  OF  THE  EQmTIONS  OF  MOTION 

The  solutions  of  the  equations  of  motion  eire  assumed  to  be  plane  waves, 
where  the  usual  cyclic  boundary  conditions  of  Bom  (25)  are  applied;  they 
insure  the  presence  of  the  characteristic  frequencies  of  the  crystal.  As  the 
body  centered  cubic  lattice  is  a  Bravais  lattice,  the  position  vector  r  of 
any  atom  in  the  lattice  may  be  expressed  as  a  linear  combination  with  integer 
coefficients  of  three  basis  vectors  a  ,  a.^,   a,  ,  This  results  in  a 
considerable  reduction  in  the  number  of  plane  wave  solutions  which  must  be 
included,  for,  if  the  motion  of  the  center  atom  is  given  by 

u  =  U  exp  27ri  (vt  +  k  •  r  ) 
o  o 

V  =  V  exp  2jt1   (vt  +  k  •  r  ) 
o  o 

w  «  W  exp  2irl   (vt  +  k  •  r  )  , 

then  the  motion  of  any  other  atom,  say  the  nth,  may  be  written  as 

u  =  U  exp  2n"i  (vt  +  k  •  r  ) 
n      *^    ^        n' 

V  =  V  exp  2ir±   (vt  +  k  •  r  ) 
n      ^    ^        n' 

w_  =  W  exp  2iri   (vt  +  k  •  r  )  . 
n  *        n 

Here  U,  V,  and  W  are  the  components  of  the  plane  wave  amplitude,  and  k  is  a 
propagation  vector  in  the  reciprocal  lattice  whose  magnitude  is  the  reciprocal 
of  the  wave  length.  The  values  of  k  are  restricted  to  the  first  Brillouin  zone, 
as  discussed  by  Brillouin  (24).  A  convenient  set  of  basis  vectors  a,,  a„,  and 
a^  is  shown  in  Fig.  k.     Using  these  vectors  one  may  write  r  as 

r   =  n^a^  +  n^i^  +  n^^j 


EXPIANATION  OF  PIATE  III 

Fig.  h.     A  set  of  basis  vectors  for  a  body  centered 
cubic  lattice. 
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where  n, ,  n^,   and  n,  are  integers. 

If,  however,  e^,  Cp,  and  e  are  unit  vectors  parallel  to  the  x,  y,  and  z 
coordinate  axis  and  a  is  the  cube  edge  length,  then  the  vector  r  may  be 
written  in  the  form  used  by  De  Launay  (22),  namely,  , 


r  e 


2  _  - 


'2x  - 


(n^  +  ^)  a  e^^  +  J-  a  e^  +  (n^  +  j-)     e^ 


(15) 


The  vector  k  may  also  be  written  with  components  k,,  kp,  and  k,  parallel  to 
the  X,  y,  z  coordinate  axis.  The  substitution  of  these  solutions  into  the 
equations  of  motion  gives  three  homogeneous  linear  equations  in  U,  V,  and  W. 
A  non-trivial  solution  requires  that  the  determinant  of  the  coefficients  of 
U,  V,  and  tf  vanish,  which  determines  the  possible  values  of  the  frequency. 
This  secular  determinant  is 


A-  +  mu) 


B 
C 


B 


Ag  +  mio 


C 
D 


A,  +  mu) 


0  , 


ilk) 


where 


A^  =  -  [8  I  +  16  7^  (^)  +  2k  72  if-)]    (1  -  cos<^  3^  cos$  ^     cos$  ^) 

2  i6 

-  4p  sin     *i  +  2  73_  (q-)   (5  cos  2   $^-cos  2  $^-cos  2  §2-cos  2  f,-2) 

rl6, 


-  3  72  (q-)  (2  +  cos  2  $^  -  cos  2  $^-cos  2  tg-cos  2  §  ) 
B  -  -  8  [|  -  7i  (~)  +  ^2  (^)  (|)3  cos  $5  sin  $^  sin  fg 
C  -  -  8  [|  -  7i  (^)  +  72  (|^)  (|)]  cos  §2  sin  «i  sin  i^ 
D  -  -  8  [2!  -  7^  (1§)  +  y^  (1§)  (|)]  eos  «^  sin  §2  sin  f^ 


(15) 


In  Eq.  (15)  the  quantity  i^   =  Trak^,  i  =  1,  2,  5,  is  used  for  convenience. 
The  quantity  a  is  the  distance  between  next-nearest  neighbors  for  the  body 


EXPIAMTION  OF  PIATE  IV 

Fig.  5.  The  first  Brillouin  zone  for  a  body  centered 

oubic  lattice  with  the  l/^J-Sth  zone  accentuated, 


PLATE  IV 
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Figure  5 
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REIATI0N3HIPS  BETWEEN  ELASTIC  CONSTANTS  AND  lATTICE  FORCE  CONSTANTS 

In  the  limit  of  long  vavelengths  where  the  wavelength  is  much  greater 
than  the  lattice  spacing  the  plane  waves  are  effectively  propagating  in  a 
continuous  medium.  This  allows  the  identification  of  the  lattice  force 
constants  with  the  elastic  constants  of  the  mediiun  found  from  measurements 
of  the  velocity  of  sound.  The  problem  then  becomes  the  formulation  of  the 
ecL^^'^ions  of  motion  governing  small  displacements  in  an  homogeneous  medium. 
This  problem  in  classical  mechanics  is  discussed  in  many  texts  on  the  subject, 
for  example,  the  necessary  equations  are  presented  by  Joos  (26), 
Theoretical  Physics,  Chapter  VIII.  The  relationship  between  strains  and 
stresses  is  given  by  Hooke's  Law 

^ij  "  J  f  ^ij-tm  ^-U  '  ^'  J'  -^^  "^  =  1'  2,  5,  (17) 

where  P..  is  the  stress  tensor,  eo   is  the  strain  tensor,  and  T. . o  is  a 
tensor  of  fourth  rank.  As  is  customEiry,  the  six  independent  components  of 
the  strain  tensor  will  be  represented  by  e.  where 

^ii  =  ^ii  '      «ij  =  2€^j         i  1^  J 

Then  from  Joos  (26),  p.  l66,  the  e       eu-e  defined  by 

_  ^u  Sv  ^W 

®11  -  3^  '     ^22  "  33F  '     ®53  "  31    ' 


u      =  ^  +  Jw        -      3w      ^u  3u  ^  Sv 

-23  "  51  +  33?  '     ^15  3^  +  3^  '    ®12  "  3^  +  35E 


> 


where  u,  v,  and  w  are  the  rectangular  components  of  the  displacement.  The 
elements  of  T  may  be  identified  with  the  elastic  constants  of  the  medium  in 
Voight's  notation.  In  this  notation  the  linear  relationship  between  stress 
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and  strain  nay  te  written  in  matrix  notation  as 

P  =  Ce 

vhere  P  is  a  row  vector  whose  elements  are  the  six  independent  elements  of  the 
stress  tensor,  and  e  is  a  row  vector  with  the  six  independent  elements  of  the 
strain  tensor  as  components.  The  matrix  C  in  general  has  21  independent 


elements  C   which  cure  related  to  the  tensor  elements  T^ 


rs  "-iila   • 

with  cubic  symmetry,  as  shown  by  Zener  (27),  C  reduces  to  the  matrix 


For  materials 


'=u. 

°12 
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0 
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0 

"iS 

=11 

^12 
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0 

"is 

=12 

^11 
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%k 
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^kk 

0 

0 

0 

0 

0 

0 

c, 

(18) 


With  the  aid  of  Eq.  (18)  the  relationships  between  the  stresses  and  the 
displacements  nay  be  found. 

Ihe  equations  of  motion  for  the  continuum  are  obtained  from  the  elastic 
equilibrium  between  the  body  forces  and  the  stresses,  which,  as  given  by 
Joos  (26),  p.  174,  are 


^  ^^^11  ,  ^^12  ^  ^ 


h?. 


12 


Sp, 


22 


SP, 


-^ 


^  "  ox 

ox      oy      oz 


g 


0  , 
0  , 
0  , 


(19) 


where  g^^,  g  ,  and  g^  are  the  components  of  the  body  forces  acting  on  a  unit 
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^2 
volume.  In  the  absence  of  external  body  forces  g  =  -p  — p  ,   the  inertial 

force  per  unit  volume,  and  upon  substitution  for  P   In  terms  of  e^^,  Eq.  (19) 


'±y 


becomes  the  wave  equation 


i.2  ;v2     n2      n2 

P  ^  -  ^^11   °12^  ^^  ^12  ^  2  +     ^     -' 


b\ 


'11  ^  -12'  ^..2  •   12  ^^^   ^^y   ^^y 


?^2    n2    >2     N2 

/OU.OU.OV   ,0W\ 

dy    ^z    3xSy   Sx^z 


(20) 


with  similar  equations  in  the  y  and  z  directions.  Assuming  plane  wave 
solutions  of  the  form 

u  =  U  exp  2iri  (vt  +  k  •  r) 
V  e  V  exp  2n'i  (vt  +  k  •  r) 
w  =  W  exp  2irl  (vt  +  k  •  r) 

results  in  the  secular  determinant 


(<3ll  -  ^kO  4  +  %k  ^^  -P  ^^   (^12  "■  %0  h^2   (°12  +  %0  V^ 

(°12  +  ^kO  h^2   (^11  -  %0  4  +  %k  ^^  -P  ^^   (^12  +  C4I.)  ^2  S 
(C^  +  Cj^j^)   kj^k^   (C^2  +  Ci^i^)  k^k^   (Oj^j^  -  C^^i^)   k|  +  Cj^i^  k^  -p  v^ 


=?  0 


(21) 


restrictiniE  the  values  of  the  frequency.  A  review  of  this  material  is  given 
by  De  Launay  (22),  p.  265. 

The  identification  of  the  lattice  force  constants  may  now  be  made  either 
from  Eq.  (20)  or  Eq.  (2l).  In  both  cases  the  limiting  condition  of  small 
relative  displacements,  i.e.  wavelengths  long  compared  to  the  lattice  spacing, 
is  made.  If,  in  Eq.  (l2)  representing  the  notion  of  the  center  atom  in  the 
x-direction,  we  expand  the  displacements  to  second  order  about  the  central 
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atom  a  coapaxison  vith  Eq.  (20)  may  be  made.  In  Eq.  (20),  for  a  body  centered 
cubic  lattice,  p  =  2m/a  ,   where  a  is  the  lattice  spacing.  An  equally  conven- 
ient comparison  may  be  made  in  the  limiting  case  of  long  wavelengths  from  the 

o 
secular  Eq.  (l^)  by  replacing  sin  $  .  by  $  .  and  cos  $  .  by  1-  $  ./2  .  In 

either  case  the  following  relationships  hold: 

a  C^^  =  2  (|  +  p  +  6  7^  +  5  7p 

a  C^2  =  2  (f  -  5  7{  +  I  7^)  (22) 


a  C^l,  =  2  (f  -.  7  •  +  I  7  •) 


'3 
where 

y{   =  (9-)  7-L   and   7^  =  (^-)  7^ 

If  7^  =  0,  for  example,  the  above  three  equations  determine  the  lattice  force 
constants  \miquely,  otherwise  the  lattice  parameters  are  determined  only  within 
a  given  ratio  of  l^l^-^'     Figure  6  shows  the  variation  of  p/a  with  ratios  of 
72/7^  for  .-the  element  vanadi\im.  The  elastic  constants  for  vanadium  are  those 
of  Alers  (28).  Solving  Eqs.  (22)  for  a,  p,  7,  in  order  to  show  the  dependence 
of  the  lattice  force  constants,  one  obtains: 

a  =  i  a  (0^2  +  5  0^4)  -  9  7^ 

P  =  ^li|-^H^i2^-B^4ii^+9  ^2  (25) 

^i  =  H  ^""hk  -  ^12^  -  ¥  ''2 

Thus  if  /g  is  zero,  7^  is  negative  for  all  elements  in  which  Cj>  is  - 
less  than  C^^  .     There  are,  however,  experimental  elastic  constants  such 
as  those  for  sodium  given  by  Quimbly  and  Siegel  (25)  at  150''K  where  C>  >  is 
greater  than  Z^^   in  which  case  7  is  positive  with  y^   equal  to  zero. 


EXPIAmTION  OF  PIATE  V 

Fig.  6.     The  variation  of  p/a  with  l-Jl^.  ^^^  ^'^^  element 
vanadivim. 
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with  the  experimental  maximvun  frequency.  Table  5  gives  the  frequency  histo- 
grams for  vanadium  for  ratios  of  yJl-^   of  0,5  and  1.0,  and  Fig.  7  shows  the 
frequency  histograms  for  these  ratios  along  with  the  case  for  /«  =  0.  Figure 
7  indicates  that  increasing  the  ratio  of  yJf\  reduces  the  power  of  the  first 
peak  and  moves  the  two  maxima  closer  together  producing  better  agreement  with 
experiment.  It  was  not  the  object  of  this  paper  to  initiate  a  parameter  study; 
however,  it  does  seem  that  the  availability  of  the  additional  force  constant 
7  may  be  used  to  alter  the  shapes  of  the  calculated  frequency  distribution  in 
order  to  obtain  better  agreement  with  experiment. 

To  better  understand  the  effect  that  the  elastic  constants  play  in  the 
determination  of  the  frequency  distribution  vanadium  distributions  were 
calculated  using  +  10^  values  for  each  elastic  constant  independently  and  with 
7^  =  0  .  The  resulting  frequency  histograms  are  given  in  Table  \,   and,  briefly, 
it  is  noted  that  increasing  C^  or  0^  by  10^  tends  to  spread  the  peaks  while 
increasing  Cj^j^  brings  the  maxima  together  slightly.  It  is  also  apparent  that 
increasing  C,,  and  O^l  by  10^  increases  the  sharpness  of  the  low  frequency  peak 
while  increasing  0      has  the  opposite  effect.  None  of  these  frequency  distri- 
butions seem  to  fit  experiment  appreciably  better  than  the  distribution  with 
the  unchanged  elastic  constants,  and  as  shown  in  Table  5,  the  calculated 
maximum  frequencies  are  all  in  the  neighborhood  of  the  experimental  maximum. 
Clark  (lO)  has  calculated  the  frequency  distribution  for  vanadium  and  the  model 
used  here  represents  the  experimental  frequency  distribution  as  well  as  the 
ordinary  central  force  model  used  by  Clark  and  predicts  a  maximiim  frequency  in 
better  agreement  with  experiment. 

In  addition  to  the  frequency  spectra  discussed  above  the  model  was  used  to 
calculate  frequency  dispersion  curves.  The  calculations  were  compared  with 
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Table  3.     Frequency  histograms   for  vanadium  with  Yo^Yi    ~  ^-^  and   1.0. 


Frequency 

'Fraction  in 

iFraction  in 

1 Frequency 

.'Fraction  in 

; Fraction  in 

interval 

! interval 

: interval 

.'intei 

rval 

: interval 

.'interval 

:y2/yi  -  0.5 

:y2/yi  =  1-0 

:y2/yi  =  0-5 

:y2/yi  =  1-0 

0. 

-  0.01 

0. 

0. 

0.50 

-  0.51 

0.02641256 

0.02621722 

0.01 

-  0.02 

0. 

0. 

0.51 

-  0.52 

0.02441675 

0.02445073 

0.02 

-  0.03 

0.00000849 

0.00000849 

0.52 

-  0.53 

0.02242944 

0.02251437 

0.03 

-  0.04 

0.00001698 

0.00001698 

0.53 

-  0.54 

0.02130839 

0.02129141 

0.04 

-  0.05 

0.00005945 

0.00005945 

0.54 

-  0.55 

0.02013639 

0.01996654 

0.05 

-  0.06 

0.00007643 

0.00007643 

0.55 

-  0.56 

0.01939752 

0.01918520 

0.06 

-  0.07 

0.00011040 

0.00011040 

0.56 

-  0.57 

0.01841235 

0.01830195 

0.07 

-  0.08 

0.00016136 

0.00016136 

0.57 

-  0.58 

0.01769046 

0.01765650 

0.08 

-  0.09 

0.00024629 

0.00024629 

0.58 

-  0.59 

0.01702803 

0.01690064 

0.09 

-  0.10 

0.00029724 

0.00029724 

0.59 

-  0.60 

0.01616177 

0.01631464 

0.10 

-  0.11 

0.00037368 

0.00038217 

0.60 

-  0.61 

0.01485388 

0.01573713 

0.11 

-  0.12 

0.00045861 

0.00045011 

0.61 

-  0.62 

0.01294300 

0.01526153 

0.12 

-  0.13 

0.00056052 

0.00055203 

0.62 

-  0.63 

0.01174552 

0.01464156 

0.13 

-  0.14 

0.00066243 

0.00067092 

0.63 

-  0.64 

0.01109158 

0.01421692 

0.14 

-  0.15 

0.00078983 

0.00078133 

0.64 

-  0.65 

0.01039517 

0.01319779 

0.15 

-  0.16 

0.00096817 

0.00096817 

0.65 

-  0.66 

0.00993656 

0.01195784 

0.16 

-  0.17 

0.00101913 

0.00098516 

0.66 

-  0.67 

0.00950342 

0.01044612 

0.17 

-  0.18 

0.00120597 

0.00123994 

0.67 

-  0.68 

0.00907029 

0.00968177 

0.18 

-  0.19 

0.00140131 

0.00132487 

0.68 

-  0.69 

0.00875606 

0.00909577 

0.19 

-  0.20 

0.00152870 

0.00156267 

0.69 

-  0.70 

0.00851826 

0.00880702 

0.20 

-  0.21 

0.00177499 

0.00173253 

0.70 

-  0.71 

0.00833991 

0.00831443 

0.21 

-  0.22 

0.00197881 

0.00197881 

0.71 

-  0.72 

0.00820403 

0.00836539 

0.22 

-  0.23 

0.00223360 

0.00219113 

0.72 

-  0.73 

0.00807664 

0.00780486 

0.23 

-  0.24 

0.00244592 

0.00237798 

0.73 

-  0.74 

0.00801718 

0,00800020 

0.24 

-  0.25 

0.00261577 

0.00267522 

0.74 

-  0.75 

0.00799171 

0.00764351 

0.25 

-  0.26 

0.00308288 

0.00298096 

0.75 

-  0.76 

0.00821252 

0,00781336 

0.26 

-  0.27 

0.00319328 

0.00318479 

0.76 

-  0.77 

0.00833142 

0,00788979 

0.27 

-  0.28 

0.00365189 

0.00350752 

0.77 

-  0.78 

0.00869661 

0,00806814 

0.28 

-  0.29 

0.00400859 

0.00388120 

0.78 

-  0.79 

0.00903632 

0.00818704 

0.29 

-  0.30 

0.00439077 

0.00428885 

0.79 

-  0.80 

0.00999600 

0,00868811 

0.30 

-  0.31 

0.00471349 

0.00466254 

0.80 

-  0.81 

0.01151621 

0.00950342 

0.31 

-  0.32 

0.00520607 

0.00509567 

0.81 

-  0.82 

0.01127841 

0.01127841 

0.32 

-  0.33 

0.00584304 

0.00566468 

0.82 

-  0.83 

0.01149923 

0,01236549 

0.33 

-  0.34 

0.00624219 

0.00610631 

0.83 

-  0.84 

0.01188140 

0.01250987 

0.34 

-  0.35 

0.00696408 

0.00664984 

0.84 

-  0.85 

0,01270521 

0.01377530 

0.35 

-  0.36 

0.00761802 

0.00731228 

0.85 

-  0.86 

0.01406405 

0.01555029 

0.36 

-  0.37 

0.00837389 

0.00795773 

0.86 

-  0.87 

0.01501524 

0.01591548 

0.37 

-  0.38 

0.00928261 

0.00897687 

0.87 

-  0.88 

0.01655243 

0.01791128 

0.38 

-  0.39 

0.01020832 

0.00967328 

0.88 

-  0.89 

0.01797073 

0,02045912 

0.39 

-  0,40 

0.01142279 

0.01072638 

0.89 

-  0.90 

0.02146976 

0,02348255 

0.40 

-  0.41 

0.01275616 

0.01177099 

0.90 

-  0.91 

0.02111306 

0,02049309 

0.41 

-  0.42 

0.01457362 

0.01338463 

0.91 

-  0.92 

0.01547385 

0,01556727 

0.42 

-  0.43 

0.01686667 

0.01539741 

0.92 

-  0.93 

0.01329120 

0.01325724 

0.43 

-  0.44 

0.01998352 

0.01745267 

0.93 

-  0.94 

0.01129540 

0.01154169 

0.44 

-  0.45 

0.02552931 

0.02175002 

0.94 

-  0.95 

0.00972424 

0.00997052 

0.45 

-  0.46 

0.03735127 

0.02511317 

0.95 

-  0.96 

0.00850128 

0.00858621 

0.46 

-  0.47 

0.03945748 

0.03113455 

0.96 

-  0.97 

0.00710846 

0.00716791 

0.47 

-  0.48 

0.02909628 

0.03848080 

0.97 

-  0.98 

0.00575811 

0.00586002 

0.48 

-  0.49 

0.02639557 

0.02910477 

0.98 

-  0.99 

0.00425489 

^m     9     ^^  ^^  ^f   %^  %^  ^^  ^^  ^^ 

0.00433132 

0.49 

-  0.50 

0.02501974 

0.02659091 

0.99 

-  1.00 

0.00219113 

0.00215717 
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dispersion  curves  for  iron  as  determined  by  Lov  (21),  Using  the  elastic  con- 
stants given  by  Low  (21)  for  iron  at  l6'C  as 

C„  =  2.352  X  lO-"-^  dyne-cm'^ 

12-2 
C, p  «  1*555  X  10   dyne-cm 

12        -2 
C. )  =  1.180  X  10   dyne-cm   , 

dispersion  curves  were  calcxilated  along  the  (lOl),  (ill),  and  (OOl)  directions. 
The  results  of  these  calculations  are  shown  in  Figs.  8,  9,  and  10.  For 
comparison  the  same  calculations  were  performed  using  a  two  constant  central 

force  model.  These  calculations  showed  that  the  shape  of  the  dispersion  curves 

12 
remained  essentially  the  same;  however,  the  maximum  frequency  of  9*755  x  10   cps 

for  the  three  constant  model  was  higher  than  the  maximum  frequency  of 

12 
8.711  X  10   cps  found  for  the  two  constant  model.  Compai-ison  with  Low's  curves 

showed  that  the  three  constant  model  represented  experiment  better  than  did  the 

two  constant  model;  however,  both  were  quite  good.  For  example  at  the  point  of 

three  fold  degeneracy  along  the  (OOl)  direction  Low  has  about  .O'^k  ev,   the  three 

constant  model  gives  .035^  ev,  and  the  two  constant  model  gives  .O325  ev.  For 

the  higher  of  the  two  transverse  modes  along  the  (lOl)  direction  at  the  zone 

boundary  the  three  constant  model  gives  .0252  ev,  the  two  constant  model 

.0232  ev,  and  Low  has  approximately  .027  ev.  At  the  point  of  three  fold 

degeneracy  along  the  (ill)  direction  the  model  employed  here  gives  .O308  ev, 

the  central  force  model  gives  .0276  ev,  and  Low  has  approximately  .O3O  ev. 

For  the  maximum  value  of  the  longitudinal  wave  along  the  (ill)  direction  Low 

has  approximately  .037  ev,  while  the  three  constant  model  gives  .O363  ev,  and 

the  two  constant  model  .O525  ev.  In  general,  the  inclusion  of  the  angular 

stiffness  terms  have  improved  the  agreement  with  experiment  in  this  case. 
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^'^  -,  ,        CONCLUSION 

The  dyimmical  properties  of  a  body  centered  cubic  lattice  of  identical 
particles  have  been  investigated  using  the  Born-von-Karman  model  of  a  crystal 
lattice.  The  model  employed  contained  non-central  as  well  as  central  forces, 
thus  the  Cauchy  relation  is  not  required.  In  particular,  the  forces  between 
lattice  points  consisted  of  non-central  angular  stiffness  forces  and  central 
forces  between  nearest  and  next -<near est  neighbors.  The  non-central  angxilar 
stiffness  force  constants  7,  and  y^  and  the  central  force  constants  Ct  and  P 
were  identified  with  the  elastic  constants  of  the  material  by  the  usual 
comparison  with  continuxun  theory. 

Frequency  distribution  histograms  have  been  constructed  for  vanadium  and 
compared  with  the  experimental  frequency  distributions  available  from  slow 
neutron  experiments.   Qualitative  agreement  was  found  in  so  far  as  the  shapes 
of  the  experimental  and  theoretical  distributions  were  concerned,  in  that  both 
distributions  showed  two  pronounced  peaks.  The  calculated  maximum  frequency 
for  vanadium  was  in  good  agreement  with  that  derived  from  the  experimental 
frequency  distributions,  which  constitutes  an  improvement  over  the  results  of 
the  central  force  model.  It  has  also  been  shown  that  the  addition  of  the 
fourth  force  constant,  7^,  could  be  used  to  alter  the  shape  of  the  frequency 
distribution  without  destroyii^g  t|ie  agreement  between  theoretical  and 
experimental  maximum  frequency.  Vanadium  distributions  were  also  calculated 
with  each  elastic  constant  varied  independently  by  +  10  per  cent.  These 
distributions  also  exhibited  qualitative  agreement  with  experiment,  and  the 
calculated  maximum  frequencies  were  generally  in  good  agreement  with 
experiment.  For  tungsten,  since  no  direct  experimental  results  are  available 


^^ 


the  calculated  dlstributioa  was  compared  with  previous  theoretical  results. 
This  calculation  was  performed  with  7^^  =  0  and  the  resulting  distribution 
conteiined  two  peaks  as  does  the  distribution  from  the  two  constant  central 
force  model;  however,  the  maximum  frequency  calculated  with  central  forces 
only  was  approximately  25  per  cent  lower  than  that  from  the  model  used  here. 
Dispersion  curves  for  iron  were  also  calculated  and  found  to  be  in  good 
qualitative  and  quantitative  agreement  with  experiment. 

The  calculations  performed  using  this  model  indicate  that  it  exhibits  a 
more  reasonable  behavior  than  the  models  generally  used  for  cubic  crystals. 
It  wo\ad  appear  that  the  improvement  is  due  in  part  to  the  removal  of  the 
restrictive  and  unrealistic  Cauchy  relation.  In  view  of  these  results  it 
would  be  of  interest  to  adapt  this  model  to  face  centered  cubic  crystals. 
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ABSTRACT 

The  Born-von  Karman  model  of  a  crystal  lattice  has  been  vised  to  Investi- 
gate the  dynamical  properties  of  a  body  centered  cubic  lattice  of  identical 
particles.  A  two-constant  model  using  central  force  interactions  between 
nearest  and  next-nearest  neighbors  has  most  frequently  been  employed  by  other 
investigators.  The  use  of  central  forces  has  the  disadvantage  of  requiring 
that  the  Cauchy  relation  be  satisfied  between  elastic  constants  of  the  crystal. 
The  model  used  in  this  work  contains  non-central  as  well  as  central  forces.  As 
a  consequence  it  has  the  inherent  advantage  of  removing  the  rfequireraent  of  the 
Cauchy  relation.  In  particular,  the  forces  between  lattice  points  consist  of 
non-central  angular  stiffness  forces  as  well  as  central  forces  between  nearest 
and  next-nearest  neighbors.  The  force  constants  employed  can  be  identified 
from  continuum  theory  in  the  usual  way  with  the  elastic  constants  of  the  lattice. 

Plane  wave  solutions  to  the  equations  of  motion  result  in  a  secular 
equation  for  the  characteristic  frequencies  of  the  lattice.  The  roots  of  this 
equation  were  found  numerically  at  approximately  59,000  points  using  an  IBM 
7090  digital  computer.  Frequency  distribution  histograms  have  been  constructed 
for  vanadium  and  compared  with  the  experimental  frequency  distributions 
available  from  slow  neutron  experiments.  Qualitative  agreement  was  found  in  so 
far  as  the  shapes  of  the  experimental  and  theoretical  distributions  were 
concerned.  The  calculated  maximum  frequency  for  vanadium  was  in  good  agreement 
with  that  derived  from  the  experimental  frequency  distributions.  This 
constitutes  an  improvement  over  the  results  of  the  simpler  central  force  model. 
Histograms  were  calculated  for  tungsten  and  compared  with  previous  theoretical 
results  since  no  direct  experimental  results  for  tungsten  are  available. 
Dispersion  curves  for  iron  were  also  calculated  and  found  to  be  in  good  agreement 


with  experiment. 

The  present  model  has  the  advantage  that  the  Cauchy  relation  is  not 
required.  In  addition,  the  calculations  performed  using  this  model  indicate 
that  it  exhibits  a  more  reasonable  behavior  than  the  model  generally  xised 
heretofore  for  cubic  crystals. 


